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Abstract: The current research proposes a novel technique 
using simulation-analytical method to determine the lifetime of 
droplet evaporation in compression ignition (Cl) engine. Using 
this technique, investigation was carried out to determine the 
optimum conditions for minimum droplet lifetime. The 
variation of droplet lifetime was studied for three parameters 
namely in cylinder pressure, compression ratio and exhaust 
gas recirculation (EGR) percentage. The best conditions for 
low droplet lifetime were obtained. 
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Introduction - 

Ignition is the initial phase of combustion. It is characterized by 
rapid increase of species temperature. The study of combustion 
and understanding of the in cylinder process enable us to design 
a much more fuel efficient engine. However, the difficulty 
encountered in combustion study is the cost involved in 
experimental setup and the complexity involved in 
understanding the mechanism. Numerical models are developed 
to understand the combustion process [1]. Commercial 
thermodynamic tools provide another approach to investigations 
in combustion. 

Effective combustion can occur only when the ignition process is 
initiated before the piston reaches TDC. Combustion process 
actually begins when the liquid fuel vaporizes completely to 
form the mixture. Once the fuel is injected into the combustion 
chamber through the injector, the liquid spray undergoes various 
physical processes and interacts dynamically with the turbulent 
fluid inside the combustor, followed by splitting up of the liquid 
fuel into droplets which collide with other droplets and coalesce 
with them. The highmomentum of the liquid spray creates a 
turbulent flow field and gas entrainment. Droplet spray may 
impinge on the cylinder wall surfaces due to the tight 
confinement inside the intake manifold or cylinders and form a 
liquid film which will eventually evaporate. Droplets can thus be 
considered the basic element which provides the fuel vapor; 
Understanding of single-droplet evaporation and combustion 
processes therefore provides important guidance in design of 
practical burners. [2] 

Droplet evaporation in the present research is considered a 
convective drop. The reason is that, when a liquid is injected in 
the combustor, there exists a relative velocity between the drop 
and the surrounding oxidizer (mostly air). It is surrounded by a 
convective boundary layer. The shape of the drop is assumed to 
be spherical. The laws of heat transfer and conservation of 
energy are then applied to this drop and analytical relationships 


are developed. Diesel-RK is used for simulation. The parameters 
used in the analytical relationship have been acquired through 
simulation. The most important parameter involved for our 
research is the expression for the lifetime of evaporation of the 
drop. The lifetime is directly related to how fast a fuel droplet 
can vaporize. Smaller droplets evaporate faster. The faster the 
droplet evaporates, the more efficient is the combustion. This in 
turn reduces the particulate matter emissions. In this research, a 
method has been developed to calculate the lifetime of an 
ethanol fuel droplet. The method developed has then been used 
to analyze the lifetime of the ethanol fuel drop under various 
design conditions with EGR. An optimum engine design 
condition is then suggested. Plots are made for the life time of 
fuel drop versus ICP and comparison is made for different CRs 
and EGR percentage. Bio-Ethanol has been chosen in this 
research as ethanol woulduse as a major fuel in the future as 
predicted by B. Hahn-Ha Gerdal et al. [4]. The properties of bio 
ethanol were obtained from ref [10]. 

II. Methodology - 

The analytical expression used to determine the life time of a 
fuel drop is derived based on the work of Me Callister et all. [2] 
as follows. The principle of conservation of energy is applied to 
the convective fuel drop. 

Energy stored in drop = energy convected. 

This can be mathematically expressed as follows, 

}hf S = JsD*.ql ( 1 ) 

Where, 

q s = heat flux convected from the droplet (W/m 2 ) 
h fg = latent heat of vaporization of droplet (J) 

D = Diameter of the fuel droplet (m) 

P/= Density of the fuel drop (kg/m 3 ) 


On solving further we get, 


dn* 

^7 = -2C&- 8 

(2) 

Where, 

1 i 


A 1 

UJcBajPraCTrt-T*) 
^ “ Pl h fg 

(3) 

And 


_ **\'T a -T b : 

(4) 


The constant Ci is assumed for simplicity purposes to be 1/2, 
i.e., the thermal layer is equal to the radius of the droplet. [2] 
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Assuming an average Reynolds number and integrating equation 
3, we have, 

D 2 = Dq — 2C ± ^ 0 t— fit (5) 

Solving the above differential equation, 

= lC^+13 (6) 

Here, is the required expression for calculating the droplet 
lifetime of evaporation. [2] 

Simulation: 

DIESEL-RK, a full cycle thermodynamic engine simulation 
software was chosen here as the results obtained from Diesel-RK 
was in close proximity to those obtained from experimental 
study. This was shown by Kuleshov A.S in ref [5,6,7]. 

The engine considered in the current research was set to the 
following parameters 


Table: 1 Experimental Engine Data 


S.No 

Parameter 

Value 

1 

Cylinder Bore diameter 

150mm 

2 

Piston stroke length 

180mm 

3 

No. of cylinders 

6 

4 

Engine Speed 

1500rpm 

5 

No. of injector 

1 

6 

Distance between spray 
centre and cylinder head 

1mm 

7 

Engine Type 

Turbocharged 


Table: 2 Variable parameters 


1 

EGR 

10%, 20%, 30%, 40% 

2 

CR 

27,28,29 

3 

Nozzle diameter(mm) 

0.195,0.200,0.210,0.225 


T able 3: Ethanol Fuel Properties 


S.No 

Parameter 

Value 

1 

Specific gravity 

0.794 

2 

Lower calorific value IMJ/kg 

27.2 

3 

Energy density (net) @ 15°C[MJ/1] 

21.6 

4 

Enthalpy of vapourisation @ b.p [kJ/kg] 

845 

5 

Research octane number 

111 

6 

Motor octane number 

89 

7 

Cetane number 

8 

8 

Stoichiometric air-fuel ratio [kg/kg] 

8.98 

9 

Spontaneous ignition temp°C 

365 

10 

Flammability limit in air (air/fuel)mass 

2.7-18.5 

11 

Flammability limit 

0.49-3.3 

12 

boiling point 

78.5 

13 

Composition 

C 2 H 5 OH 


The ICP vs 0, T vs 0 , injection velocity versus 0, mass fraction 
of fuel vaporization vs crank angle and outlet fuel velocity 
versus crank angle plot were obtained from the simulation. The 
numerical values pertaining to the plots were obtained using 
Microsoft Excel 2007. The T and ICP values were noted at a 0 
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of 13° before TDC. The simulation also gives the sauter diameter 
of droplet (d32). The plots obtained for one case is shown in Fig 
1,2, 3, 4, 5. The mixture formation and combustion simulation is 
also shown. 



Fig 1. ICP (bar) versus O ( degree ) (CR 28 with 10%EGR) 



Fig 2. T (K) versus Q (degree) (CR 28, 10%EGR) 



Crank Angle after Injection Beginning (degree) 

Fig 3. Injection Velocity (m/s) versus G (degree) (CR 28, 
10%EGR) 
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Figure :4 Variation fuel velocity outlet of Nozzle with Crank 
Angle ( rpm 1500, cr 28: 1) 



Figure :5 Variation of Evaporated fuel mass fraction with 
Crank Angle ( rpm 1500, cr28: 1) 
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FIG: 6 Mixture formation and Combustion Simulation - 
Compression ratio 28:1 


Table :4 Combustion Simulation Data( cr28:l) 


No 

Parameter 

Expansion 

Value 

1 

A/F_eq 

Air Fuel Equivalence Ratio in 
the Cylinder 

2.7070 

2 

F/A_eq 

Fuel Air Equivalence Ratio in 
the Cylinder 

0.36941 

3 

p_max 

Maximum Cylinder Pressure, 
bar 

217.84 

4 

T_max 

Maximum Cylinder 
Temperature, K 

1562.4 

5 

CA_p.max 

Angle of Max. Cylinder 
Pressure, deg. A.TDC 

2.0000 

6 

CA_t.max 

Angle of Max. Cylinder 
Temperature, deg. A.TDC 

15.000 

7 

dp/dTheta 

Max. Rate of Pressure Rise, 
bar/deg 

11.510 

8 

p_inj .max 

Max. Injection Pres, (before 
nozzles), bar 

1603.6 

9 

d_32 

Sauter Mean Diameter of 
Drops, microns 

8.9534 

10 

Theta_i 

Injection / Ignition Timing, 
deg. B.TDC 

13.000 

11 

Phi_inj 

Duration of Injection, deg 

26.500 

12 

x_e.id 

Fuel Mass Fraction 
Evaporated during Ignit. 
Delay 

0.05879 

13 

Rs_tdc 

Swirl Ratio in the Combustion 
Chamber at TDC 

0.10000 

14 

Rs_ivc 

Swirl Ratio in the Cylinder at 
IVC 

0.06194 


Estimating the life time of drop: 

The parameters found using the simulation method are then 
substituted in the following relationships. Average Temperature 
(T avg ): The ambient temperature is assumed to be 303K. The 
average temperature is hence determined. 

~ 2 


[12] Average thermal conductivity of air is given by 

k = Cl. 5207 0- u ):r 3 - C4.S574ff- Ije 3r : +- 
Cl .0184s -m )T — 3 ,9333 s -B4 

The viscosity of air is given by [8] 

^ ^ lr 0 / T + S 

Where, pu = (1.716)10 _s — , T 0 =273.1K, =110.56 K. 

m-s 
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I CP 

Air density is calculated by [9]p = —Where, p is the 

density (kg/m 3 ), ICP is the In cylinder Pressure (Pa), R(J/Kg-K) 
is the gas constant and Tis the In cylinder temperature(K). 

M-Cp 

Prandtl Number is given by [3]Pr = — 

Specific heat capacity [11] given by: 

C p = (1.9327 9~ w )T 4 - (7:9'9'9'9e-° 7 )r a -f 
(1,1407 9~ m )T 2 - (4 .48 9 0 e “ DL ) 7 T 1.0575 9 m 

III. Results and Discussion- 

The droplet life time corresponding to the various input 
parameters are tabulated. Graphs have been plotted to show the 
variation of the droplet lifetime with ICP for various CRs and 
EGR percentage. 


Ta ble: 5 Droplet Life Time Variation with EGR 10% 


CR 27 

CR 28 

CR 29 

ICP 

f [L/j 

ICP 

hifa 00 

ICP 

hifa 00 

(bar) 

(x 10-*) 

(bar) 

(x 10 - *) 

(bar) 

10-*) 

94.74 

714 

98.15 

686 

101.37 

662 

94.92 

506 

98.25 

487 

101.43 

470 

94.96 

345 

98.28 

332 

101.46 

320 

94.97 

317 

98.37 

305 

101.48 

295 


Ta ble: 6 Droplet Life Time Variation with EGR 20% 


CR 27 

CR 28 

CR 29 

ICP 

(bar) 

00 
(X 10- 4 ) 

ICP 

(bar) 

*Bfe 00 
(x 10- 4 ) 

ICP 

(bar) 

t[L/a 00 

(X 10-*) 

96.53 

693 

99.78 

668 

103.09 

644 

96.55 

492 

99.84 

474 

103.11 

458 

96.57 

335 

99.86 

323 

103.19 

312 

96.6 

308 

99.93 

297 

103.23 

287 


ISSN:2319-6890)(online), 2347-50 13 (print) 

1 August 2016 

The above values are plotted in a graph. They are shown as 
follows: 



Fig 7 : Droplet Lifetime Vs. ICP (EGR 10%) 



In Cylinder Pressure 

Fig :8 Droplet Lifetime Vs ICP (EGR 20%) 


Table: 7 Droplet Life Time Variation with EGR 30% 


CR 27 

CR 28 

CR 29 

ICP 


ICP 

^lifa G 5 ) 

ICP 

hife 00 

(bar) 

(x 10-*) 

(bar) 

(X 10- 4 ) 

(bar) 

(x 10-*) 

93.19 

694 

96.19 

669 

99.62 

645 

93.28 

492 

96.52 

475 

99.67 

458 

93.3 

335 

96.54 

323 

9968 

312 

93.31 

308 

96.57 

297 

99.71 

286 


Ta ble: 8 Droplet Life Time Variation with EGR 40% 


CR 27 

CR 28 

CR 29 

ICP 

Hfa 00 

ICP 

hifz 00 

ICP 

hifz 00 

(bar) 

(x 10-*) 

(bar) 

(x 10-*) 

(bar) 

(x 10-* 

89.08 

700 

92.24 

675 

95.27 

651 

89.09 

497 

92.26 

479 

95.27 

462 

89.14 

338 

92.26 

326 

95.275 

315 

89.19 

311 

92.27 

300 

95.28 

290 



Fig: 9 Droplet Lifetime Vs ICP (EGR 30%) 
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In Cylinder Pressure 


Fig: 10 Droplet Lifetime Vs ICP (EGR 40%) 

It can be seen from the graphs that the droplet lifetime decreases 
with compression ratio and incylinder pressure. From Table (5) 
for an EGR of 10% for a compression ratio of 27, the value of 
doplet lifetime has decreased from 714e-04 to 317 e-04 s. This 
indicates a drop of about 55.6% in the droplet lifetime. Now 
keeping the ICP constant and observing the change in 
compression ratio, it is observed that the droplet life time drops 
by about 6.9% from a CR 27 to CR 29. For an EGR of 20%, the 
droplet life time has dropped by 55.43% by increasing ICP for 
CR29 and a 6.81% drop by increasing the CR from 27 to 29 at 
constant ICP. For an EGR of 30%, there is a decrease of 
55.503% in droplet lifetime for CR 29 and decreasing ICP, and 
when ICP is constant, there is a drop of 0.681% by increasing 
CR from 27 to 29. With 30%EGR, there is a drop of 55.45% in 
droplet lifetime for CR29 and increasing ICP, while there is a 
6.752% drop when ICP is fixed but CR is increased from 27 to 
29. Hence, it can be inferred that as the compression ratio 
increases and as the in cylinder pressure which indirectly refers 
to fuel injection pressure increases, the droplet life time 
decreases. This is because, with higher injection pressure, the 
fuel is injected at a higher velocity as the pressure energy at the 
injector nozzle is converted to kinetic energy. Now, when the 
fuel is injected at such high pressures, the fuel spray will split 
into droplets of very small diameters; a smaller drop will 
naturally evaporate faster as shown in fig 7- 16. A low 
evaporation time indicates a better mixing of the fuel with the 
oxidizer in the combustion chamber. This results in a higher 
proportion of the fuel being oxidized giving a better combustion 
efficiency. Also, this indicates a lower emission of particulate 
matter. 

IV. Conclusion: 

The method developed to determine the droplet lifetime was 
developed from fundamental principles [3] and also validates the 
claim that combustion characteristics improve by increasing 
compression ratio [13] and fuel injection pressure [14]. Also, a 
Cl engine with EGR 30% and CR29 and nozzle injector diameter 
0.222 mm gives the lowest droplet lifetime of vaporization. 


Nomenclature: 

T 

ICP 

EGR 

CR 

T-Life 

0 


In Cylinder Temperature 
In Cylinder Pressure 
Exhaust Gas Recirculation 
Compression Ratio 
Droplet Lifetime 
Crank Angle 
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